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Abstract 

Background: Fibromyalgia is associated with substantial socioeconomic loss and, despite considerable research 
including numerous randomized controlled trials (RCTs) and systematic reviews, there exists uncertainty regarding 
what treatments are effective. No review has evaluated all interventional studies for fibromyalgia, which limits 
attempts to make inferences regarding the relative effectiveness of treatments. 

Methods/design: We will conduct a network meta-analysis of all RCTs evaluating therapies for fibromyalgia to 
determine which therapies show evidence of effectiveness, and the relative effectiveness of these treatments. We 
will acquire eligible studies through a systematic search of CINAHL, EMBASE, MEDLINE, AMED, HealthSTAR, 
PsychlNFO, PapersFirst, ProceedingsFirst, and the Cochrane Central Registry of Controlled Trials. Eligible studies will 
randomly allocate patients presenting with fibromyalgia or a related condition to an intervention or a control. 
Teams of reviewers will, independently and in duplicate, screen titles and abstracts and complete full text reviews 
to determine eligibility, and subsequently perform data abstraction and assess risk of bias of eligible trials. We will 
conduct meta-analyses to establish the effect of all reported therapies on patient-important outcomes when 
possible. To assess relative effects of treatments, we will construct a random effects model within the Bayesian 
framework using Markov chain Monte Carlo methods. 

Discussion: Our review will be the first to evaluate all treatments for fibromyalgia, provide relative effectiveness of 
treatments, and prioritize patient-important outcomes with a focus on functional gains. Our review will facilitate 
evidence-based management of patients with fibromyalgia, identify key areas for future research, and provide a 
framework for conducting large systematic reviews involving indirect comparisons. 

Keywords: Fibromyalgia, Systematic review, Network meta-analysis, Multiple treatment comparison, Randomized 
controlled trials 



Background 

Fibromyalgia is a syndrome characterized by chronic 
widespread pain and excessive tenderness at 11 of 18 spe- 
cific muscle-tendon sites, for which no clear cause can be 
found [1]. Approximately 2% of the general population in 
the United States suffers from fibromyalgia, with women 
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affected ten times more often than men [2]. Similar preva- 
lence rates have been reported in Canada (33%), Brazil 
(4.4%) and Western European countries, including 
Germany (3.2%), Spain (2.4%), Italy (2.2%), Sweden (2.5%), 
France (1.4%), Italy (3.7%) and Portugal (3.6%) [3-8]. 
Treatment directed towards fibromyalgia is highly variable 
and long-term prospective observational studies have 
found that patient outcomes, even in specialized rheuma- 
tology clinics, are typically poor [9,10]. 

Poor treatment outcomes have led to frustration 
among patients and their clinicians, and contributed to 
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the impact of fibromyalgia on disability [11-13]. Obser- 
vational studies have found that 20% to 50% of persons 
with fibromyalgia report that they are unable to work or 
that they can work only a few days per month, and 27% 
to 55% receive disability or social security payments. Of 
those that do work, 36% experience an average of two or 
more absences from work per month [11-13]. 

Physicians typically find fibromyalgia difficult to man- 
age and patients with fibromyalgia are likely to be dissat- 
isfied with treatment [14]. A recent survey of 1,200 
primary care physicians in the United States (33% re- 
sponse rate) found that only 14% of respondents indi- 
cated very good or excellent satisfaction with managing 
patients with fibromyalgia and other medically unex- 
plained symptoms [15], and a study of 400 British gen- 
eral practitioners (75% response rate) found that only 
44% of respondents felt there were effective treatment 
options available for this population [16]. 

There have been no less than 38 systematic reviews 
addressing therapies for fibromyalgia [17-54]; however, 
the large majority of reviews to date have explored ther- 
apies in isolation or looked at a subgroup of treatments. 
There have been two network meta-analyses we are 
aware of that explored multiple therapies for fibromyal- 
gia [51,52]. The first [52] restricted comparisons to 
pharmacological treatments only and used the Jadad 
scale to assess study quality [55], which has been criti- 
cized as overly simplistic and placing too much em- 
phasis on blinding [56]. The second [51] excluded 
complementary and alternative medicine approaches, 
had no assessment of study quality, and summarized 
continuous data using the standardized mean difference, 
which is vulnerable to differential variability in popula- 
tions enrolled and adds challenges in interpreting the 
magnitude and importance of treatment effects [57,58]. 
No review has looked at all interventional studies for 
fibromyalgia, which limits attempts to make inferences 
regarding the relative effectiveness of available treat- 
ments. Further, none used the GRADE framework to 
establish confidence in pooled estimates of treatment ef- 
fect [59]. 

We will explore all therapies for fibromyalgia that have 
been tested in randomized controlled trials (RCTs) and 
use Bayesian mixed treatment comparison methods 
(adjusted indirect comparisons) to complement the dir- 
ect comparisons of the relative effects of competing 
interventions in RCTs [60,61]. 

Methods/design 

Protocol and registration 

Our protocol is registered on PROSPERO (CRD4201200 
3291), http://www.crd.york.ac.uk/PROSPERO. 

Our paper conforms to the PRISMA guidelines for 
reporting systematic reviews [62]. 



Eligibility criteria 

Trials eligible for our review will (1) have enrolled adult 
patients (>18 years of age) presenting with fibromyalgia 
or a related condition (for example, myofascial pain 
syndrome, fibrositis, fibromyositis, muscular rheuma- 
tism, chronic generalized pain syndrome), and (2) have 
randomized patients to an intervention or a control arm. 

Information sources and search 

We will identify relevant RCTs, in any language, by a sys- 
tematic search of CINAHL, EMBASE, MEDLINE, AMED, 
HealthSTAR, PsycINFO, PapersFirst, ProceedingsFirst, 
and the Cochrane Central Registry of Controlled Trials, 
from inception of the database, with relevant MeSH 
headings. An experienced librarian (QM) developed a 
sensitive search strategy for each individual database 
(see Additional file 1). We will scan the bibliographies 
of all retrieved trials and other relevant publications, 
including reviews and meta-analyses, for additional 
relevant articles. 

Study selection 

Twenty reviewers with experience in health research 
methodology will work in pairs to screen, independently 
and in duplicate, titles and available abstracts of identi- 
fied citations and acquire the full text publication of any 
article that either reviewer judge as potentially eligible. 
The same reviewer teams will independently apply eligi- 
bility criteria to the full text of potentially eligible trials. 
Reviewers will resolve disagreements by consensus or, if 
a discrepancy remains, through discussion with one of 
two arbitrators (JWB or GHG). 

Data collection process and data items 

Using standardized forms (see Additional file 2) and a 
detailed instruction manual that will be used to inform 
specific tailoring of an online data abstraction program 
(DistillerSR), ten teams of reviewers will extract data in- 
dependently and in duplicate from each eligible study. 
To ensure consistency across reviewers, we will con- 
duct calibration exercises before starting the review. 
Data abstracted will include demographic information, 
methodology, intervention details, and all reported 
patient-important outcomes. Reviewers will resolve dis- 
agreements by discussion, and one of two arbitrators 
(JWB or GHG) will adjudicate unresolved disagree- 
ments. We will contact study authors to resolve any 
uncertainties. 

Two reviewers (GC, EAC) will independently extract 
details on interventions and outcomes from all RCTs in 
order to classify them into common intervention categor- 
ies. The abstractors will develop the categories independ- 
ently and then achieve consensus through discussion. 
Outcomes will be classified into domains, based on the 
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guidelines published by the Initiative on Methods, Meas- 
urement, and Pain Assessment in Clinical Trials 
(IMMPACT) [63-70]. These consist of nine domains: pain; 
physical and emotional functioning (including Quality of 
Life); participant rating of improvement and satisfaction 
with treatment; adverse events; participant disposition (for 
example, adherence to the treatment regime and reasons 
for premature withdrawal from the trial); role functioning 
(that is, work and educational activities, social and recre- 
ational activities, home and family care); interpersonal 
functioning (that is, interpersonal relationships, sexual 
activities); sleep; and fatigue. 

Risk of bias in individual studies 

Reviewers will assess risk of bias within each study with a 
modified Cochrane risk of bias instrument which assesses 
the following key domains: sequence generation; alloca- 
tion concealment; blinding of participants, healthcare 
professionals, outcome assessors, data collectors, and 
data analysts; incomplete outcome data; selective out- 
come reporting; and other sources of bias. Reviewers 
will input response options of 'definitely yes', probably 
yes', probably no', and 'definitely no' for each of the do- 
mains, with 'definitely yes' and probably yes' ultimately 
assigned low risk of bias and 'definitely no' and 'prob- 
ably no' assigned high risk of bias [71]. Reviewers will 
resolve disagreements by discussion, and one of two 
arbitrators (JWB or GHG) will adjudicate unresolved 
disagreements. 

Direct comparisons meta-analyses 

To pool outcome data for trials that compare the same 
intervention with the same comparator, we will use 
random effects meta-analyses, which are conservative in 
that they consider both within and among study differ- 
ences in calculating the error term used in the analysis 
[11]. We will pool cross-over trials with parallel design 
RCTs, using methods outlined in the Cochrane hand- 
book to derive effect estimates [72]. Specifically, we will 
perform a paired £-test for each cross-over trial if either 
of the following are available: 1) the individual partici- 
pant data; 2) the mean and SD (or standard error) of the 
participant-specific differences between the intervention 
and control measurements; 3) the mean difference and 
one of the following: (i) a t-statistic from a paired £-test, 
(ii) a p-value from a paired £-test, or (iii) a confidence 
interval from a paired analysis; or 4) a graph of measure- 
ments of the intervention arm and control arm from 
which we can extract individual data values (pending 
that the matched measurements for each individual can 
be identified) [72]. If these data are not available, we will 
approximate paired analyses by first calculating the 
mean difference for the paired analysis (MD = M E - M c ) 
and the standard error of the mean difference: 



SE(MD) = — where N represents the number of par- 
ticipants in the trial, and SD^ represents the standard de- 
viation of within-participant differences between the 
intervention and control measurements [72]. If the stand- 
ard error or standard deviation of within-participant differ- 
ences is not available, we will impute the standard deviation 
using methods outlined in the Cochrane Handbook [73]. 

Ensuring interpretable results from pooled estimates of 
effect 

We will use a number of approaches to facilitate inter- 
pretable results from our meta-analyses. For trials that 
report dichotomous outcomes, we will calculate the odds 
ratio (OR) to inform relative effectiveness. We will 
acquire estimates of baseline risk from observational 
studies or, if not available, from the median of the con- 
trol group from eligible RCTs. 

When pooling across trials that report continuous out- 
comes using the same instrument, we will calculate the 
weighted mean difference (WMD), which maintains the 
original unit of measurement and represents the average 
difference between groups [74]. Once the WMD has 
been calculated, we will contextualize this value by 
noting the corresponding minimally important differ- 
ence (MID) - the smallest change in instrument score 
that patients perceive is important. We will prioritize 
use of anchor-based MIDs when available, and calculate 
distribution-based MIDs when they are not. 

Contextualizing the WMD through the MID can be 
misleading because clinicians may interpret a WMD less 
than the MID as suggesting that no patient obtains an 
important benefit, which is not accurate. Therefore, we 
will generate an estimate of the proportion of patients 
who have benefited by applying the MID to individual 
studies, estimating the proportions who benefit in each 
study, and then aggregate the results in order to provide 
a summary estimate of the proportion of patients who 
benefit from treatment across all studies. Details of the 
methods by which we will conduct this analysis are 
presented immediately below in our discussion of situa- 
tions in which investigators have used different instru- 
ments to measure the same construct. 

For trials that use different continuous outcome mea- 
sures that address the same underlying construct, one 
cannot calculate a weighted mean difference, and we will 
therefore calculate a measure of effect called the 
standardized mean difference (SMD) or 'effect size'. This 
involves dividing the difference between the intervention 
and control means in each trial (that is, the mean differ- 
ence) by the estimated between-person standard devi- 
ation (SD) for that trial. The SMD expresses the 
intervention effect in SD units rather than the original 
units of measurement, with the value of a SMD 
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depending on both the size of the effect (the difference 
between means) and the SD of the outcomes (the inher- 
ent variability among participants). 

This common approach to pooling continuous out- 
come data is often problematic. If the heterogeneity of 
patients is different across studies, the SD will vary 
across studies. Therefore, given the same true difference 
in outcome between intervention and control groups, 
trials with more heterogeneous patients will show appar- 
ently - but spuriously - smaller effects than trials enrolling 
less heterogeneous patients. Furthermore, interpretation 
of the magnitude of effect when represented as SD units is 
not intuitive. 

In order to address these issues, we will contextualize 
the SMD value through MID units, which are not vul- 
nerable to the distortions that varying heterogeneity of 
populations can create and are more interpretable to 
both clinicians and patients [58,75]. For outcome mea- 
sures that have an established anchor-based MID we will 
use this measure to convert the summary effect into OR. 
We will complement this presentation by either 
converting the summary effect into natural units of a 
widely accepted instrument used to measure changes in 
the domain of interest or, if such an instrument is not 
available, we will substitute the MID for the SD (denom- 
inator) in the SMD equation, which will result in more 
readily interpretable MID units instead of SD units [58]. 
Finally, we will, as for SMD, provide a summary estimate 
of the proportion of patients who benefit from treatment 
across all studies 

We illustrate this approach with the following 
example. We will first describe how we will summarize 
the outcome in MID units. Assume that a trial reports a 
mean difference (MD) on a continuous outcome meas- 
ure "X", and assume that an anchor-based MID for 
instrument X, MID X , has been established. The esti- 
mated MD is a random variable. If we standardize this 
random variable by dividing it by the MID X , we get a 
new random variable, MD/MID X . We know from basic 
probability theory that because MID X is a constant, the 
variance of MD/MID X is given by: 



Var 



MD \ _ Var[MD) 



MID\) MID\ 



That is, the variance of the mean difference divided by 
the square of the MID. Further, the standard error of 
MD/MID X is given by: 



SE 



MD 
M1D\ 



IVar(MD) _ SE(SD) 
MID\ MID X 



Consider a meta-analysis that included k trials. The 
first j trials use disease-specific instrument A, and the 



last k-j trials use disease-specific instrument B. Let MD t 
denote the mean difference observed in trial i, let MID A 
denote the MID established for instrument A, and let 
MID B denote the MID established for instrument B. Fur- 
ther, let mi denote the MID standardized effect for trial i. 
To pool results across trials using MIDs we must first 
estimate the m t and its associated variance for all trials. 
For i = 1, . . ., j we have: 



MDt VariMDi) 
and Varym-i) - 



MID A 

and for i = j + 1, . . ., k we have: 



MDj 
MID B 



and Var (mi) 



MIDI 



Var(MDi) 
MID\ 



By defining the trial weights as Wi = Var(mi)" 1 , we can 
use the fixed-effect model inverse variance method to 
pool the MID-standardized mean differences using the 
formula: 



m 



Where m denotes the pooled MID-standardized mean 
difference. The standard error of m can be calculated 
using the formula: 



se(m) = 1/(5^ 



Wi 



The associated confidence intervals can subsequently 
be derived. MID-standardized mean differences can also 
be combined in a random-effects model using weights 
Wi = (Var(mi) + t 2 )" 1 , where t 2 denotes the between-trial 
variance. 

This presentation does not address the risk that clini- 
cians may interpret all mean effects below the MID as 
unimportant, and presume important benefit for all pa- 
tients when mean effects exceeds the MID. We will ad- 
dress this issue by assuming normal distributions of data 
and then calculating the proportions of participants in 
the intervention and control groups in each study that 
demonstrated an improvement greater than the MID 
[76]. The results are then pooled across studies. If we 
only have post-test data (rather than magnitude of 
change), we will apply this approach if evidence exists 
regarding meaningful thresholds. For instance, if one 
knows that people with scores of less than 8 on the 
Hamilton rating scale for depression (HAM-D) are con- 
sidered to be not depressed, one could examine the pro- 
portion of individuals below that threshold. 

If such meaningful thresholds do not exist, we will use 
post-test data and assume that the minimally important 
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change within an individual corresponds, on average, to 
the minimally important difference between individuals. 
Making this assumption, one can calculate the difference 
in the proportion who benefit in intervention and control 
To do this, we will take the mean value in the control 
group plus one MID unit, and calculate the proportion of 
patients in each group above that threshold. 

If an anchor-based MID has not been established for 
all instruments, we will assume a meta-analysis control 
group probability (p c ) and use the SMD to calculate the 
OR. Specifically, we will construct a conceptual meta- 
analysis control group with mean standard deviation 
a c > and group size n c , and a conceptual meta-analysis 
intervention group with mean p E , standard deviation o~ E , 
and group size n E such that the SMD can be represented 
as SMD = p E - [i c and o E = o c = 1. We will set /u c = 0 and 
(?c = l> and our threshold (T) will be equal to 0' 2 (pcX 
where O' 1 is the inverse standard normal cumulative 
distribution function. We will then use the derived 
threshold to calculate the conceptual intervention group 
probability (p E ). The intervention group mean response 
is assumed to follow a normal distribution with mean 
SMD and a SD of 1. Thus, the intervention group prob- 
ability is p E = 1 - 0(T- SMD), Having estimated the 
conceptual meta-analysis control and intervention group 
probabilities from the pooled SMD, we will calculate the 
OR as follows: 

0R= Mi-pc) 

PC{1 -PE) 

To calculate the 95% CI, we will use the above formu- 
las substituting the upper and lower bounds of the SMD 
confidence interval. We will complement this presenta- 
tion by converting the SMD into natural units of a 
widely accepted instrument used to measure changes in 
the domain of interest or, if such an instrument is not 
available, we will calculate the ratio of means [58]. 

Assessment of heterogeneity and subgroup analyses 

We will examine heterogeneity of meta-analyses using 
both a chi-squared test and the I 2 statistic, the percent- 
age of among-study variability that is due to true differ- 
ences between studies (heterogeneity) rather than 
sampling error (chance) [77,78]. We have generated the 
following a priori hypotheses to explain variability be- 
tween studies: (1) interventions will show larger effects 
among trials enrolling >50% patients clearly defined as 
meeting the American College of Rheumatology (ACR) 
criteria for fibromyalgia [1] versus trials in which less 
than 50% of subjects meet the ACR criteria; (2) larger ef- 
fects in trials of patients in which receipt of disability 
benefits or involvement in litigation was an exclusion 
criteria versus those that did not exclude patients on this 



basis; and (3) studies with greater risk of bias will have lar- 
ger effects than studies with lower risk of bias. We will 
perform subgroup analyses on a component-by-compo- 
nent basis if we detect variability within the individual risk 
of bias components. We will conduct z-tests [79] to estab- 
lish if subgroups differ significantly from each other [80]. 

Multiple comparison meta-analysis 

We will examine the assumptions of similarity (for an 
indirect treatment comparison) and consistency (for 
mixed-treatment comparison) before conducting a net- 
work meta-analysis. To assess relative effects of compet- 
ing treatments, we will construct a random effects 
model within the Bayesian framework using Markov 
chain Monte Carlo methods in WinBUGS (MRC Bio- 
statistics Unit, Cambridge, UK) [81]. We will model out- 
comes in every treatment group of every study and 
specify the relations among the effect sizes across studies 
[82]. This method combines direct and indirect evidence 
for any given pair of treatments. We will use the 
resulting 95% credible intervals (Crls) to assess treat- 
ment effects [83]. We will assess the probability that 
each treatment was the most efficacious regimen, the 
second best, the third best, etcetera, by calculating the 
effect size for each treatment compared with an arbitrary 
common control group and counting the proportion of 
iterations of the Markov chain in which each treatment 
has the highest effect size, the second highest, etcetera. 
A key assumption behind multiple treatments meta- 
analysis is that the analyzed network is consistent; that 
is, that direct and indirect evidence on the same com- 
parisons do not disagree beyond chance. We will locate 
and estimate inconsistencies by employing a mixed 
treatment comparisons inconsistency model in the 
Bayesian framework [84]. 

We will use a recent users guide published in the 
Journal of American Medical Association (of which two 
of our team were authors: KT and GHG) to assess the 
strength of inferences and credibility of our network 
meta-analysis [85]. We will use this guide to critically 
appraise our work in 3 domains and 12 subdomains 
(Table 1). 

Confidence in pooled estimates of effect 

Reviewers will, independently and in duplicate, assess 
the confidence in effect estimates for all outcomes using 
the GRADE (Grading of Recommendations, Assessment, 
Development and Evaluation) rating system [86]. In the 
GRADE system of rating quality of evidence for each 
outcome, randomized trials begin as high quality evi- 
dence, but may be rated down by one or more of five 
categories of limitations: (1) risk of bias, (2) consistency, 
(3) directness, (4) imprecision, and (5) reporting bias 
[80]. After considering these categories, the confidence 
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Table 1 Critical appraisal guide for a Network Meta-Analysis [85] 

Critical appraisal criteria Specific items 

A. Are the results of the study valid? 

Did the review explicitly address a sensible clinical question? 
Was the search for relevant studies exhaustive? 
Were there major biases in the primary studies? 

B. What are the results? 

What was the amount of evidence in the network? 

Were the results similar from study to study? 

Were the results consistent in direct and indirect comparisons? 

What were the overall treatment effects and their uncertainty, and how did the treatments rank? 
Were the results robust to sensitivity assumptions and potential biases? 

C. How can I apply the results to patient care? 

Were all patient-important outcomes considered? 

Were all potential treatment options considered? 

Are any postulated subgroup effects credible? 

What is the overall quality and what are limitations of the evidence? 



in estimates for each outcome will be categorized as 
follows: 'high' quality of evidence (we are very confident 
that the true effect lies close to that of the estimate of 
the effect); moderate' quality of evidence (we are moder- 
ately confident in the effect estimate and the true effect 
is likely to be close to the estimate of the effect, but 
there is a possibility that it is substantially different); 
low' quality of evidence (our confidence in the effect 
estimate is limited and the true effect may be substan- 
tially different from the estimate of the effect); and Very 
low' quality of evidence (we have very little confidence 
in the effect estimate and the true effect is likely to be 
substantially different from the estimate of effect) [86]. 

We will assess publication bias by visually observing 
asymmetry of the funnel plot for each outcome. As a 
rule of thumb [87], tests for funnel plot asymmetry 
should be used only when there are at least ten studies 
included in the meta-analysis. Otherwise the power of 
the tests is too low to distinguish chance from real 
asymmetry. We will report our results by type of inter- 
vention (for example, psychotherapy, analgesics, antide- 
pressants) and focus on patient-important outcomes. 
We will report all direct comparison data, and only 
complement these data with indirect comparison data if 
the strength of inferences from direct comparisons is 
similar or less than the strength of inferences from indir- 
ect comparisons. If there are no direct comparisons 
available, we will report indirect comparison data only. 

Knowledge translation 

We plan to create a stakeholder advisory committee with 
representation from ambulatory health care providers 
from across Ontario, Canada as well as from key organiza- 
tions. We will ensure that we have geographically diverse 



representation including primary care providers who prac- 
tice in rural areas of the province. Members of our stake- 
holder committee will be invited to attend our planning 
meeting and share their input/advice with members of the 
review team. 

Our team also will engage in an end-of-study knowledge 
translation workshop. The purpose of this activity will be 
to share our findings with key relevant stakeholders 
(researchers, clinicians and decision-makers) in order to 
identify future opportunities for dissemination, beyond 
traditional peer-reviewed publications, with our stake- 
holders, discuss how to maximize uptake of our findings in 
patient education and clinical practice, and determine fu- 
ture research directions. The overall goal of the workshop 
is to develop an agenda that will establish directions to 
develop and implement our research findings into practice. 

The following strategies will be used to promote 
awareness of the stakeholder meeting findings according 
to the Ottawa Model of Research Use in which informa- 
tion is tailored to specific audiences: (1) distribution of 
findings to all involved participants for further input, 
sharing within their organization, and for sharing with 
their own stakeholders via newsletter, web site, or other 
methods; (2) presentation at relevant peer-reviewed and 
community conferences; and (3) publication in an open- 
source peer- reviewed journal. We anticipate that this 
meeting will identify new areas of inquiry for research 
and practice, such as the development of new educa- 
tional tools for patients and clinicians. We also antici- 
pate that new collaborations and networks will be 
created that will support the identified work going for- 
ward. Any groups identified through the meeting will be 
included as part of the report back to the stakeholders 
in order to broadly disseminate the findings. 
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Discussion 

Our review will evaluate all treatments for fibromyalgia, 
provide relative effectiveness of treatments, and evaluate 
the quality of the evidence in a thorough and consistent 
manner using the GRADE approach [88-90]. Addition- 
ally, many existing reviews focus on surrogate outcomes, 
such as number of tender points, stiffness, range of mo- 
tion, or laboratory values; we will prioritize patient- 
important outcomes such as function and quality of life. 
The results of our systematic review will be of interest 
to a broad audience including patients diagnosed with 
fibromyalgia, health professionals managing fibromyal- 
gia, employers, human resource professionals, insurers/ 
compensation boards, and labor groups. Our review will 
facilitate evidence-based management of patients with 
fibromyalgia, identify key areas for future research, and 
provide a framework for conducting large systematic re- 
views involving indirect comparisons. 

Additional files 



Additional file 1: Appendix A. Search Strategy. 
Additional file 2: Appendix B. Data Abstraction Forms. 
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